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T  h  ACT 


n  limited  area  sever,  layer  physical -numerical  model  Tor  the  Lower 
trcpospiieric  rerior.  (surface-lbOOm)  is  described.  The  grid  interval 
is  naif  triat  of  tire  standard  numerical  weather  prediction  grid  used  in 
the  hemispneric,  free  atmospheric,  operational  model  at  the  Air  iorce 
llobal  Weather  Central  (Alois).  Tnis  model  is  an  integral  part  of  the 
complete  nil..!,  ueso-scale  (sub-synoptic)  numerical  analysis  and  pre¬ 
diction  system.  This  model  provides  greater  horizontal  and  vertical 
resolution  in  both  the  numerical  analyses  and  numerical  forecasts.  It 
is  used  to  predict  tne  more  detailed  smaller  scale  atmospheric  per¬ 
turbations  which  are  important  in  specifying  sensible  weather  elements. 


Important  features  of  this  boundary  layer  model  include:  a 
pletelv  automated  objective  numerical  analysis  of  input  data;  the 


a  com- 


trans:. ort  of  heat  and  moisture  by  three  dimensional  wind  flow  (including 
terrain  and  fractionally  induced  vertical  motions);  latent  heat  ex- 
cuange  in  water  substance  phase  changes:  and  eddy  flux  of  heat  and 
water  vapor. 


xnput  data  are  conventional  synoptic  surface  and  upper  air  re¬ 
ports.  other  prediction  models  provide  horizontal  wind  components 
at  the  upper  boundary  and  an  estimate  of  cloudiness  above  the  boundary' 
layer  region,  forecasts  for  the  lower  boundary7  and  surface  layer  are 
empirically  derived,  despite  some  approximations  which  broadly  simplify 
the  real  planetary  boundary  layer  processes,  operational  use  indicates 
the  model  is  capable  of  producing  detailed  forecasts  out  to  24  hours. 


winter  case  study  is  discussed. 
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This  memorandum  describes  the  boundary  Layer  Model  being  used  by 
the  Air  force  Global  Weather  Central  (AFGv.'C)  to  produce  automated  fore¬ 
casts  of  temperature,  specific  humidity  and  tliree  dimensional  wind 
flow  in  the  lower  1600  meters  of  the  troposphere.  This  model  has  been 
in  operational  production  since  March,  1969. 

The  material  is  presented  in  five  sections.  First,  some  background 
information  is  presented  to  place  the  model  in  perspective  as  to  its 
development  and  integration  into  the  A:  GV.'C  system,  next,  a  mathematical 
and  physical  description  of  the  model  is  given.  Third,  the  procedures 
used  to  obtain  initial  and  forecast  boundary  conditions  are  presented, 
fourth,  a  description  of  the  computational  model  is  presented,  and, 
filially,  a  case  study  is  shown. 

There  were  several  important  factors  which  acted  as  guiding  or 
limiting  influences  on  the  development  of  this  model.  The  need  for  an 
operational  model  which  would  treat  the  thermo-hydrodyamic  processes 
in  the  planetary'  boundary  layer  has  existed  for  many'  years.  Any  model 
which  would  handle  the  strong  horizontal  and  vertical  gradients  over 
highly  variable  terrain  with  appropriate  consideration  of  adiabatic 
and  diabatic  influences  as  well  as  moisture  sources  and  sinks  must 
necessarily'  be  very  complex;  thus,  computer  size  and  speed  became  im¬ 
portant  in  determining  the  mathematical  sophistication  of  the  boundary 
layer  model  as  well  as  the  size  of  the  forecast  region.  The  availability 
of  adequate  data  samplings  required  to  form  initial  conditior.  s  for  the 


m 


trendy  important  influence  upon  the  mode  l  1  s  development  (especially 


upon  the  choice  of  grid  size) 


window"  concept  using  a  fine  mesh  grid  has  beer,  in  existence 


at  AbilAJ  for  several  years.  This  svstem  uses  a  detailed  numerical 


analysis  over  a  subregion  or  limited  area  of  high  operational  interest, 


The  Ailri.  1'  Mine -Mesh  Upper  Air  Analysis  Model  (reference  i)  uses  a  grid 


interval  of  approximately  100  nautical  miles,  or  one-half  the  standard 


111  Northern  hemispheric  grid  spacing.  This  analysis  increases  the 


resolution  to  a  point  where  initial  identification  can  be  made  of 


those  smaller  scale  perturbations  so  important  to  the  forecasting  of 


sensible  weather  elements.  The  AF1WC  Keso-Scale  Prediction  Model  (re¬ 


ference  2)  is  used  to  provide  free  atmospheric  forecasts  of  up  to  2k 


hours  over  this  "window  area" . 


There  are  numerous  theoretical  boundary  layer  models  which  have 


been  used  experiment ally  to  treat  boundary  layer  phenomena,  i.e.,  sea 


breeze,  slojie  winds,  etc.  These  models  are  usually  restricted  by  the 


availability  of  input  data,  computer  running  time,  or  their  inability 


to  converge  to  a  stable  solution.  The  AFGv/C  i  oundary  Layer  Model  is 


the  first  known  operational  model  to  systematically  treat  many  thermo- 


hydrodynamical  interactions  occurring  in  the  atmospheric  boundary  layer. 


Many  approximations  had  to  be  made  to  achieve  an  operational,  model. 


Gome  of  these  approximations  are  very  simple,  but  the  ability  of  this 


model  to  simulate  and  predict,  with  reasonable  accuracy,  the  behavior 


of  the  atmosphere  in  the  boundary  layer  region  appears  to  justify 


tneir  use .  Discussions  of  those  approximations  will  be  given  as  various 


aspects  of  the  model  are  described. 
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The  A:  l.'J  boundary  Layer  Model  has  evolved  frcm  Gerrlty’s  (re¬ 
ferences  c  and  9)  numerical  Model  for  the  Prediction  of  Synoptic- 
scale  how  Cloudiness,  wliich  was  developed  under  U3AF  contract.  The 
emphasis  was  changed  from  prediction  of  low  cloudiness  to  the  predic¬ 
tion  of  the  basic  meteorological  parameters  of  temperature,  moisture 
and  three  dmnensional  wind  flow,  initial  analyses  and  forecasts  of 
these  parameters  provide  a  forecaster  with  the  basic  ingredients  not 
only  for  forecasts  of  trends  of  low  cloudiness,  but  also  for  the  prog¬ 
nosis  of  triggering  mechanisms  for  severe  weather,  low  level  mechanical 
turbulence,  inversion  layers  (for  the  specification  of  refractive  index), 
haze  layers  and  many  other  derived  meteorological  parameters. 

The  model  is  designed  and  programmed  to  use  the  AFGWC  Fine-Mesh 
analysis  as  the  first  guess.  The  model  uses  the  wind  forecasts  of 
the  AKrt.C  Meso-Scale  Prediction  Model  for  upper  boundary  conditions. 

An  important  feature  of  the  model  design  is  the  allowance  of  systematic 
investigation  of  the  boundary  layer  processes  and  continued  development. 
More  sophisticated  techniques  will  be  introduced  into  the  model  as 
more  in  learned  atout  the  real  atmosphere  within  the  boundary  layer 
and  those  physical  processes  which  couple  the  micrometeorological  and 
synoptic  scales. 


2 .  DESCRIPTION  OF  THE  MODEL 

The  AFGh'C  boundary  Layer  Model  is  a  hybrid  which  uses  diagnostic 
techniques  to  obtain  fields  of  motion  and  forecast  techniques  to  solve 
the  tendency  equations  for  temperature  and.  moisture  .  The  momentum 
equations  were  not  included  in  the  model  because  the  large  horizontal 
and  vertical  gradients  and  the  smaller  scale  perturbations  found  near 
the  ear'-:  1  s  surface  are  inadequately  described  by  observed  data.  The 
increased  computer  core  and  computing  time  necessary  to  solve  the  mo¬ 
mentum  eciuations  were  not  warranted  without  some  assurance  that  tills 
procedure  would  provide  better  wind  forecasts  than  could  be  obtained 
through  the  diagnostic  procedures.  The  diagnostic  procedures  of  this 
model  make  optimum  use  of  the  output  of  proven  free  atmospheric  numer¬ 
ical  models  to  provide  horizontal  wind  flow  at  the  upper  boundary.  The 
horizontal  wind  at  each  lower  level  is  obtained  by  using  the  detailed 
thermal  structure  of  the  boundary  layer,  geophysical  parameters  valid 
at  fixed  geographical  locations  and  computed  frictional  forces  required 
to  obtain  a  resultant  balanced  flow. 

The  mathematical  formulation  of  the  model  will  be  presented  bj- 
giving  the  canonical  form  of  the  tendency  equations  for  temperature  and 
moisture,  and  showing  how  values  are  obtained  for  each  terr..  These 
equations  are  not  unique;  derivations  may  be  found  in  references  6  and 
11.  The  tendency  equations  for  temperature,  specific  humidity  and 
specific  moisture  are : 
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where 


T  -  temperature 

u  --  horizontal  wind  component  in  the  positive  x-direction 

u  ~  horizontal  wind  component  in  the  positive  y-direction 

-  frictionally  induced  vertical  motion 
it-  --  terrain  induced  vertical  notion 

h  -  specific  humidity  (vapor) 

R  -  specific  moisture  (vapor  plus  liquid  moisture) 


R  coefficient  of  eddy  conductivity 


coefficient  of  eddy  diffusivity 


gravitational  acceleration 


f  i  ■  specific  heat  at  constant  pressure 


—  differential  of  specific  humidity  due  to  change  of 

state 


the  horizontal  wind  components  ( [(  f  )  at  each  level  are  obtained 
from  the  solution  of  the  following  modified  Kkman  equations: 
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■where  a  =  geostroprac  wind  component  in  the  positive  x-direction 

d 

v  =  geostrophic  wind  component  in  the  positive  y-direction 

y 

;■  =  height  of  level  for  which  wind  component  is  valid 

1  -  height  of  the  second  level  (50m) 

LI  -  wind  component  in  the  positive  x-direction  computed 
from  surface  layer  relations 

-  wind  component  in  the  positive  y-direction  computed 
from  surface  layer  relations 

f  n  i'/z 

a 

l  2  Kw  ■ 

•'  -  coriolis  parameter 

k  -  coefficient  of  eddy  viscosity 

Frictionally  induced  vertical  motion  ( tc )  i~  obtained  (neglecting 
variations  of  density^  by  integrating  the  continuity  equation  in  the 
form: 


.k  . v-'. 


Trie  terrain  induced  vertical  motion  (to)  is  obtained  by  solving 
the  following  equation  for  each  level: 

H  >  E 

to  »  u  —  +  u  —  (2.7) 

U  ')  IJ 

where  E  =  terrain  elevation. 

Specific  humidity  (0)  is  calculated  from  dew  point  temperature 
by  using  Teten's  equation  to  obtain  vapor  pressure  and  by  using  a 
standard  atmospheric  profile  for  the  applicable  ambient  pressure. 
Teten's  equation  for  the  relation  between  saturation  vapor  pressure 
and  temperature  is  given  by 

=  6.11  x  l oaT / (T+b) 

where  T  is  the  temperature,  in  degrees  Celsius,  a  and  b  are  con¬ 
stants.  Over  water,  a  =7.5  and  b  ~  237-3- 

Specific  moisture  { R )  is  assumed  to  equal  specific  humidity  at 
the  initial  time.  This  restriction  will  be  relaxed  after  tests  are 
run  on  the  new  AFGWC  multi-layered  nephanalyses  and  suitable  relations 
between  clouds  and  liquid  moisture  content  have  been  established. 

The  exchange  coefficients  for  temperature  and  moisture  for  each 
layer  are  computed  from  equations  developed  by  Estoque  (reference  6). 
These  relate  the  exchange  coefficients  to  stability,  wind  shear  and 


~iixing  length. 


A  determination  is  made  of  whether  free  or  forced  con¬ 


vection  exists,  .ej -ending  upon  the  i-'ici.ardson  number.  The  equation 
for  forced  convection  [R{  -  0 .  3  )is 


S'  1 


k  =  (/»;■«<) 


(2.8) 


wnere 


i.  barman's  constant  (0.38) 
empirical  constant  (-2.0) 


wind  soeed  shear  through  the  layer 


i  ichardson  number,  whicii  is  defined  by 


f  .1 


]  I  .Si 


(2.9) 


where  I  =  is  the  mean  temperature  for  the  layer 
,  =  is  the  dr}r  adiabatic  lapse  rate 

l'he  exchange  coefficient  for  free  convection  (R(  -0.3)  can  be 

written  as 


KU  =  Kll 


a  i  'i 


(2.10) 


wnere 


--  is  an  empirical  const, ant  (0.9). 


16 

The  magnitudes  of  the  eddy  coefficients  are  bounded  between  10  -  10 

cm^  sec  .  it  is  anticipated  that  these  bounds  will  be  changed  to  a 
lower  range  after  more  testing  has  been  accomplished.  The  value  of  the 
eddy  visosity  (Km)  is  computed  for  the  surface  layer  and  is  assumed 
independent  of  height  for  solution  of  equations  (2.4)  and  (2.5).  The 
exchange  coefficients  for  heat  and  moisture  are  assumed  equal  and  are 
evaluated  for  eacn  layer  of  the  model  at  each  time  step. 


(V-' 


The  horizontal  wind  components  ( U  >  t' )  at  50m  are  computed  from 
equations  applicable  only  to  the  shallow  layer  adjacent  to  the  earth's 
surface  (the  contact  layer").  An  excellent  derivation  and  discussion 
of  these  equations  can  be  found  in  reference  9.  The  calculated  values 
°1  are  functions  of  coriolis  force,  surface  roughness,  geo- 

strophic  wind  and  empirical  results  from  Blackadar  (reference  4)  which 
determine  a  frictional  velocity  and  turning  angle  for  the  contact  layer. 

The  geostrophic  wind  components  [u  ;t  vl?)are  obtained  at  each 
level  by  solving  the  equations 
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lower  bound.uri  conditions  of  J  emperature  (  T  )  raid  specific  humidity 
(P  ,j-.*  i.ecueu  be  solve  file  i  enuenc;.  equations  [2.  j_j,  {2.2)  and  (2.p). 

The  change  in  surface  temperature  at  instrument  shelter  height  is  deter¬ 
mined  by  the  use  of  climatologicaily  derived  diurnal  surface  temperature 
change  curves  for  clear  skies  and  light  winds.  These  diurnal  curves 
are  also  used  to  simulate  radiational  cooling  within  the  lower  levels 
of  the  boundary  layer.  The  temperature  change  curves  are  mean  monthly 
values  for  each  grid  point  (a  curve  for  each  month  at  each  grid  point). 
This  change  is  then  modified  by  the  amount  of  cloudiness  within  and 
above  the  boundary  layer.  Advection  and  adiabatic  changes  are  used  to 
modify  the  diurnal  temperature  change.  Finally,  surface  condensation 
is  allowed,  which  results  in  modification  of  temperature  through  the 

i 

release  of  latent  heat.  Surface  specific  humidity  forecasts  are  made 
using  advective  and  water  substance  phase  changes  only. 

The  saturation  adjustments  for  temperature,  specific  humidity  and 
specific  moisture  stem  directly  from  the  work  of  McDonald  (reference  10). 
The  technique  allows  condensation  to  occur  with  super  saturation,  and 
evaporation  to  occur  when  liquid  moisture  is  present  and  the  air  is 
not  saturated.  Precipitation,  either  through  or  from  the  boundary  layer, 
is  not  considered  in  the  present  model. 

The  coordinate  system  employed  is  a  quasi-Cartesian  system  derived 
for  a  spherically  shaped  earth  (small  terms  have  been  neglected).  This 
particular  Boundary  Layer  Model,  which  is  currently  used  for  the  North 
American  window,  has  a  29  x  27  grid.  This  grid  is  a  subset  of  the  North 
American  fine-mesh  window  and  has  a  grid  spacing  of  approximately  100 


nautical  miles  (one-half  that  of  the  standard  AFGWC  grid).  A  horizontal 
depiction  of  the  forecast  grid  area  is  shown  in  Figure  1. 
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The  inner  rectangle  encloses  the  forecast  grid,  and  the  outer  rectangle 
encloses  the  boundary  layer  analysis  grid.  If  we  consider  the  grid  at 
80  degrees  west  longitude,  the  x-coordinate  is  positive  to  the  east 
and  the  y-coordinate  is  positive  to  the  north.  The  model  includes 
eight  levels  as  shown  in  Figure  2.  The  heights  of  the  eight  levels 
are:  surface,  30,  150,  300,  600,  900,  1200,  and  1600m  above  the  surface. 
The  vertical  coordinate  j -  )  is  the  height  above  the  terrain.  The  fixed 
geometrical  thickness  between  any  two  levels  insures  that  a  detailed 
description  of  the  boundary  layer  can  be  obtained  ever  variable  terrain. 
The  size  of  the  grid,  the  number  of  levels  and  the  thicknesses  between 
levers  can  be  changed  easily  in  the  operational  model. 

Output  data  for  each  level  are  stored  hourly,  from  zero  to  24  hours 
in  the  AFGv.'d  data  base  (magnetic  drums).  The  parameters  stored  are  the 
three  components  of  the  wind,  temperature,  D-value,  geometric  height, 
specific  humidity  and  relative  humidity. 

D-values  for  each  level  are  computed  using  output  from  the  free 
atmospheric  prediction  model  to  obtain  hourly  forecasts  of  D-values 
at  1600m  above  the  terrain.  Values  at  lower  levels  are  computed  based 
on  the  specific  temperature  anomalies  between  layers. 

The  probability  that  the  boundary  layer  model  will  accurately  pre¬ 
dict  the  future  state  of  the  atmosphere  is  directly  related  to  the 
accuracy  and  representativeness  of  the  initial  and  forecast  boundary 
conditions.  The  next  section  describes  in  detail  the  techniques  used 
in  obtaining  initial  boundary  conditions  for  the  model  equations. 


3  .  ICn.'uAnY  LA  YE!-:  AI.ALY51S 

The  difficulty  of  obtaining  initial  distributions  of  temperature 
and  moisture  which  are  geometrically  and  dynamically  consistent  is 
greater  for  the  boundary  layer  region  than  for  the  free  atmosphere, 
iactcrs  contributing  to  tiiis  difficulty  are:  the  boundary  layer  region 
is  characterized  by  strong  vertical  and  horizontal  gradients;  the  finite 
difference  technique  used  to  solve  partial  differential  equations  re¬ 
quires  data  at  fixed  grid  intervals,  but  observations  are  taken  at 
arbitrary  locations;  surface  observations  may  be  for  a  height  above  or 
below  the  "model'  terrain;  the  number  of  surface  observations  over  the 
North  American  forecast  grid  is  approximately  ten  times  greater  than 
the  number  of  upper  air  observations;  and  most  observations  taken  in 
mountainous  regions  are  from  valley  reporting  stations. 

An  automated  technique  was  designed  and  programmed  to  provide  the 
best  possible  initial  numerical  analyses  of  temperature  and  specific 
humidity  with  due  consideration  of  the  difficulties  noted  above. 

The  basic  rationale  of  the  analysis  technique  is  to  analyze  from 
data  dense  regions  toward  data  sparse  regions.  In  the  case  of  the  boundary 
layer  model,  the  data-dense  level  is  at  the  surface  and  the  data-sparse 
levels  are  in  the  region  above  the  surface.  Only  conventionaa  data, 
i.e.,  surface  and  upper  air  reports,  are  used  in  the  analysis.  These 
analyses  are  based  on  a  modified  barnes  technique  (reference  2). 

The  analysis  technique  is  data  oriented;  i.e.,  it  passes  through 
the  data  list  to  determine  the  data  density  from  which  an  appropriate 
scan  radius  for  each  data  point  can  be  calculated.  This  radius  varies 
from  tvjo  to  seven  grid  distances  depending  on  data  density.  If  there 


are  no  data  within  seven  grid  distances,  the  first  guess  field  value 
for  the  grxd  point  in  question  is  used  in  the  analysis. 

The  boundary  layer  analysis  programs  are  designed  to  make  optimum 
use  of  data  while  maintaining  vortical  and  lateral  consistency.  The 
analysis  proceeds  from  the  surface  level  upward  to  the  highest  level 
cf  the  boundary'  layer  model.  This  is  accomplished  by  first  obtaining 
the  best  possible  automated  surface  analysis  and  building  upward  layer 
by  layer.  Analysis  of  lapse  rates  between  levels  rather  than  analysis 
of  actual  values  of  the  parameter  in  the  build-up  procedure  insures 
vertical  consistency  throughout  the  depth  of  the  boundary  layer. 

First  guess  fields  are  obtained  from  the  output  of  the  AFGWC 
Fine-Mesh  Analysis  Model  (reference  l).  This  output  includes:  sea  level 
pressure,  surface  temperature  and  dew  point  temperature  fields;  temp¬ 
erature  and  specific  humidity  fields  (calculated  from  dew  point  temp¬ 
eratures)  at  850,  700,  500  and  300  mb.  These  fields  are  used  to 
calculate  first  guess  lapse  rates  for  each  layer  of  the  model.  A  de¬ 
tailed  radiosonde  (RA0B)  analysis  for  each  station  over  the  analysis 
region  is  performed  to  obtain  data  corresponding  to  geometric  heights 
Ox  each  level  of  the  model.  Lapse  rates  of  temperature  and  specific 
humidity  for  each  layer  are  computed  from  each  RAOB.  This  provides 
first  guess  fields  and  observational  data  for  use  in  the  analyses. 


The  remainder  of  this  section  will  describe  how  the  data  are 


used. 


Surface  data  reported  in  RA0BS  are  incorporated  into  the  first 
guess  surface  analysis  field  in  the  following  manner.  Deviation  data 
(the  difference  between  RAOB  data  and  the  first  guess  analysis  at  the 
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a  at  a  point  location),  are  first  computed.  ...e  v^.^e  •  f  fur 


analysis  at  the  data,  location  is  found  bj  quadratic  interpolation. 


acn  deviation  influences  a  nu.  lie r  of  gr  id  points  in  the  vicinity  of 


the  data.  The  amount  of  influence  or  weight  varies  with  the  distance 


between  the  grid  point  and  the  data  point,  and  the  scan  radius  computed 


for  that  data  point.  The  weight  ( (p '  given  a  deviation  applicable  at 
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where  R  -  distaiice  from  data  point  to  grid  ’point 


SR  -  scan  radius  computed  from  data  density 


Tiie  individual  contribution  of  each  piece  of  deviation  data  at  a  grid 


pc.  Lnt  is  the  product  of  the  deviation  and  weight  computed  for  the  data 


point . 


vny  grid  point  may  be  influenced  by  several  pieces  of  data.  Thus, 


the  total  change  applied  at  a  grid  point  is  made  up  of  normalized  values 


of  the  individual  changes.  Each  piece  of  deviation  data  contributes 


to  tiie  total  change  by  an  amount  which  is  equal  to  the  product  of  the 


deviation  and  the  ratio  of  the  weight  for  each  individual  data  point 


to  the  sum  total  weight  for  all  deviations  having  influence  at  the  grid 


point.  This  method  of  weighting  the  deviation  data  helps  to  prevent 


the  formation  of  unrealistically  sharp  gradients  between  sparse  and 


dense  data  areas  (such  as  sea  coasts). 


After  the  contributions  of  ail  deviations  have  been  considered 


to  form  a  change  field,  the  first  guess  field  and  the  change  field 


are  added.  The  resulting  field  is  an  imp'roved  estimate  of  the  analysis 
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field. 


-'de  deviation  data  are  again  computed  from  the  refined  analysis 


ana  the  original  At.;  data.  The  sajne  procedures  are  followed  to  compute 


a  new  change  field.  r.'he  procedure  is  repeated  three  times,  with  each 


successive  scan  providing  closer  agreement  between  data  and  analysis. 


The  analyses  of  temperature  and  specific  humidity  above  the  sur¬ 


face  are  produced  in  a  somewhat  different  maimer  than  at  the  surface. 


i’iie  important  differences  lie  in  the  method  of  obtaining  the  first 


guess  fields  for  the  deviation  analyses  and  in  the  use  of  lapse  rate 


analyses  instead  of  actual  values  of  the  respective  parameters. 


The  first-guess  lapse  rate  analyses,  computed  from  the  output  of 


the  :■  Ine-Kesh  Upper  Air  Analysis  Ilodel,  tend  to  be  noisy,  especially 


over  variable  terrain.  Therefore,  a  certain  amount  of  smoothing  is 


necessary.  A  smooth  first  guess  field  is  obtained  by  talcing  the  actual 


1.A0:  lapse  rate  data,  spreading  the  data  to  grid  points,  and  combining 


the  resultant  field  with  the  preliminary  guess  field.  The  actual  RAOt 


lapse  rate  data  are  used  to  compute  deviation  data,  which,  in  turn, 


are  used  to  compute  change  fields  in  the  same  manner  as  at  the  surface. 


The  first  guess  field  is  then  refined  by  addition  of  the  change  field. 


This  procedure  is  repeated  three  times. 


The  magnitude  of  the  temperature  lapse  rate  is  limited  by  a  value 


cc:;  puted  from  the  equation 


1  +  L 


(3.2) 


where 


—  ' 

cut-off  value  (deg  ta  “) 


dry  adiabatic  lapse  rate  (0.0098  deg  m 


M 


i 


m 


Ml* 


.  .  -  auto convective  lapse  rate  (0.0312  deg  m  ) 

ft 

L  =  a  parameter  varying  from  or.e  to  I.,  where  one  is  the 
bottom  layer  and  F  is  the  top  layer 
'.a  thickness  of  the  layer  (m) 

The  cuo-off  value  rapidly  approaches  the  dry  adiabatic  lapse  rate  as 
the  height  above  the  surface  increases.  For  temperature  lap3e  rates, 
the  analysis  at  each  grid  point  is  checked  against  the  cut-off  value 
after  each  scan.  If  the  lapse  rate  exceeds  the  cut-off  value,  the 
lapse  rate  is  replaced  with  the  limiting  value. 

Analyses  of  temperature  and  specific  humidity  are  produced  for 
the  level  at  the  top  of  each  layer  by  adding  the  value  of  the  lapse 
rate  through  the  layer  to  the  value  of  the  parameter  at  the  next  lower 
level.  In  the  specific  humidity  analyses,  a  check  is  made  to  avoid 
negative  and  supersaturated  values.  The  analyzed  values  obtained  for 
specific  humidity  are  used  as  initial  values  of  specific  moisture . 
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/+.  THE  COMPUTATIONAL  MODEL 


xhe  previous  sections  in  this  report  dealt  mainly  with  equa¬ 
tions  used  in  the  model  and  the  method  of  analysis  used  to  provide 
initial  and  forecast  boundary  conditions.  This  section  describes  the 
design  of  the  computer  program  and  the  basic  techniques  used  to  solve 
the  equations  of  the  model. 

A  packed  data  format  is  used  to  reduce  computer  core  requirements. 
The  packing  scheme  and  data  ranges  are  given  in  Table  1.  This  packing 
scheme  allows  the  use  of  79  fieids  in  core  at  one  time  (the  same  amount 
of  the  core  would  accomodate  only  22  fields  of  unpacked  data). 

The  tendency  equations  (2.1),  (2.2)  and  (2.3)  we re  transformed  to 
a  linearized  set  of  finite  difference  equations  in  which  the  coefficients 
are  assumed  constant,  for  a.  given  time  increment.  These  equations  are 
solved  by  an  uncentered  (upstream)  difference  technique.  A  combination 
of  a  forward  difference  for  the  time  derivative  and  upstream  difference 
for  the  3pace  derivative  is  used;  this  procedure  ha3  the  advantage  of 
remaining  computationally  3tabie  as  long  as  the  standard  numerical 
stability  criterion  . T  •  is  satisfied.  Here,  .  T  is  the  time 

step  (30  minutes),  /  x  is  the  grid  interval  and  C  is  the  speed  of 
the  advecting  wind  plus  the  rate  of  movement  of  the  system.  From  a 
theoretical  viewpoint,  the  upwind  difference  approximation  is  less 
accurate  than  a  centered  difference  approximation.  However,  from  a 
practical  viewpoint  the  upwind  method  is  in  better  agreement  with  the 
physical  concept  of  advection.  in  determining  the  time  rate  of  change 
due  to  advection  of  a  quantity  at  a  given  point  in  space  and  time,  only 
the  upstream  values  f  .  qu.-a.city  re  .  f  in;  urbane*. 
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Certain  boundary  conditions  are  applied  in  the  solution  of  the 
tendency  equations,  cut  flow  is  allowed  at  the  .Lateral  boundaries  of 
the  forecast  grid  iut  inflow  is  not.  this  restriction  causes  a  de¬ 
gradation  of  the  forecast  on  the  windward  side  of  the  boundary  layer 
forecast  region,  this  degradation  penetrates  the  region  with  the 
speed  of  the  wind.  Various  tecnniques  have  been  investigated  to  re¬ 
move  this  restriction,  but  no  real  breakthrough  is  anticipated  until 
lateral  boundary  conditions  can  be  forecast  with  higher  fidelity. 

Tiie  general  faow  of  tne  complete  model  is  depicted  in  Figure  3- 
The  large  quantity  and  variety  of  data  required  to  prepare  the  initial 
and  forecast  boundary  conditions,  and  the  complexity  of  the  forecast 
model  require  sequential  operation  of  eight  main  computer  programs. 

There  are  numerous  subprograms  associated  with  each  of  these  programs. 
Six  programs  are  required  to  pre-process  and  analyze  the  data  prior 
to  execution  of  the  forecast  program,  tne  program  is  used  for  the 
prognostic  model  and  one  program  is  used  to  accomplish  the  post-  proces¬ 
sing  and  storage  of  selected  parameters  into  the  AFGUC  drum  data  base. 

The  remainder  of  this  section  will  describe  the  functions  performed 
by  each  program.  The  purpose  of  this  description  is  to  clarify  what 
data  arc  used  and  how  the  model  operates.  The  AFGUC  computer  system 
is  drum  oriented;  temporary,  immediate  access,  drum  storage  is  used. 

The  data  used  by  the  model  are  aval...  able  in  the  AFGWC  drum  data  base. 

4.1  ouund ary  layer  Analysis  Programs 

The  first  analysis  program  cnecks  the  AbGhC  dram  data  base  for  the 
presence  of  tne  data  required  to  execute  the  model.  If  the  required 
fields  a T'5  available ,  a  terr  ain  elevation  field  is  read  from  the  data 
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base.  The  field  applicable  to  the  boundary  layer  window  is  extracted 
and  written  on  a  drum  file  for  later  use  in  the  analysis  programs. 

The  second  analysis  program  reads  the  fine-mesh  analysis  fields 
from  the  drum  data  base  and  computes  fields  of  specific  humidity  for 
the  surface,  £50,  700,  500,  and  300  mb  levels.  Surface  pressure,  temp¬ 
erature  and  specific  humidity  fields  for  the  window  are  written  on 
drum  for  later  use.  finally,  this  program  computes  temperature  and 
specific  humidity  lapse  rates  for  each  layer  of  the  model.  These 
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first  guess  fields  are  also  written  on  a  drum  for  later  use  in  the 
analyses . 

The  third  analysis  program  reads  the  KACi  data  for  the  boundary 
layer  window  from  the  AFGUC  drum  data  base.  A  hydrostatic  check  is 
made  and  geopotentiau  heights  for  each  significant  level  are  computed 
using  the  hypsometric  equation.  The  reported  values  of  temperature 
and  computed  values  of  specific  humidity  for  each  sounding  are  then 
used  to  interpolate  values  for  each  model  level.  Finally,  these 
values  are  written  on  drum  for  later  use. 


The  fourth  analysis  program  reads  the  first  guess  fields  and 
KACfc  data  from  the  drum  files  and  produces  analyses  of  temperature  and 
specific  humidity  for  each  level  of  the  model. 

The  fifth  analysis  program  reads  the  analysis  fields,  fixed  field 


data  (i.e.  radiational  coefficients,  surface  roughness  and  terrain 
height)  and  cloud  data  from  drum  files,  and  arranges  the  data  in  a 
form  that  can  be  used  by  the  forecast  program.  These  data  are  packed 
into  the  appropriate  arrays.  These  arrays  are  written  on  drum  and  are 
ultimately  read  and  used  by  the  forecast  model. 


V 


The  forecast  model  program  reads  initial  conditions  from  drum 


files,  solves  the  finite  difference  versions  of  the  diagnostic  equations 


(2.11)  and  (2.12)  and  computes  the  geostropliic  wind  at  each  level.  The 


finite  difference  versions  of  equations  (2.4)  and  (2.5)  are  solved  to 


obtain  a  balanced  horizontal  flow  for  each  level.  The  finite  difference 


versions  of  equations  (2.6)  and  (2.7)  are  solved  to  obtain  vertical 


velocities  induced  by  surface  friction  and  terrain.  At  this  point. 


the  model  begins  its  prognostic  cycle  by  forecasting  the  surface  temperature 


and  surface  specific  humidity  for  the  next  time  step.  The  finite  dif¬ 


ference  versions  of  the  tendency  equations  (2.1),  (2.2)  and  (2.3)  are  then 


solved  to  obtain  temperature,  specific  humidity  and  specific  moisture 


for  the  next  forecast  time.  The  latent  heat  computation  and  appropriate 


adjustments  of  temperature,  specific  humidity  and  specific  moisture 


are  made.  The  program  is  then  ready  to  read  in  a  new  upper  boundary 


wind,  make  a  diagnosis  of  the  flow  at  the  remaining  levels  from  the 


new  forecast  temperature  distribution,  make  new  surface  forecasts,  etc, 


f 


delected  lore cast  -.iis|  .Lays  arc;  output  on  a  higli  speed  printer  for  de¬ 


sired  forecast  Intervals. 


4.3  I  ost-tlocesslng  and  storage  Programs 


The  model  out: uts  the  data  listed  in  Table  1,  at  hourly  forecast 


interv-  ,  a  drum  file.  These  data  are  packed  in  the  arrays  named 


in  the  tab^e,  using  the  horizontal  coordinates,  L  and  M,  and  the  vertical 


coordinate,  as  subscripts  identifiying  the  grid  point  and  level. 


The  last  ;  rogran  reads  and  unpacks  these  arrays,  calculates  1600m 


Ao'rL  b-values  from  the  output  of  the  AHIV/C  Meso-Scale  Prediction  Model 


preference  2),  and  subsequently  computes  relative  humidities  and  D-values 


for  each  level  of  the  boundary  layer  model.  The  data  are  then  stored 


in  the  nUniC  drum  data  base . 


Tiie  next  section  is  a  study  of  a  representative  case,  showing  some 


of  the  model's  capabilities  and  samples  of  the  available  output. 
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Tiiis  section  is  devoted  to  a  case  study  for  i6  January  1970. 


The  purpose  of  this  study  is  to  demonstrate  the  detailed  initial 


analysis,  the  applicability  of  the  diagnostic  wind  procedures  for  flow 


over  variable  terrain  at  both  initial  and  forecast  times  and  to  show 


the  type  of  output  which  is  routinely  available  from  the  model.  The 


pe.1  formance  of  tne  model  on  tills  particular  case  isn't  any  better  or 


worse  than  can  be  expected  on  a  particular  day. 


The  general  synoptic  situation  at  U000Z  on  16. January  1970  is 


depicted  in  Figure  4  (reproductions  of  facsimile  charts  transmitted 


tne  national  Meteorological  center).  The  main  storm  track  was 


through  the  Southern  United  States  as  evidenced  by  the  occluded  system 


approacliing  the  West  Coast,  the  remains  of  an  occluded  system  along  the 


hew  Mexico-Texas  border  and  the  wave  disturbance  off  Southeastern 


Florida.  Arctic  and  modified  arctic  air  dominated  Canada  and  the 


horthern  United  States.  The  cold  arctic  high  over  Western  Canada  had 


been  nearly'  stationary  during  the  previous  24  hours. 


Figures  5  thru  8  show  the  rapid  wave  development  and  formation 


of  an  occlusion  in  the  Great  Lakes  Region,  the  cold  penetration  of 


arctic  air  along  the  Eastern  slopes  of  the  Rocky  Mountains,  the  develop¬ 


ment  of  a  wave  along  the  arctic  front  in  Eastern  Colorado  and  the  move¬ 


ment  of  the  occluded  system  into  the  Western  United  States. 


a  few  streamlines  have  bee.,  placed  oi  f  igure  9  to  show  the  wind  flow 


at  loOOm  above  the  terrain.  This  wind  field  was  derived  from  neignt 


fields  from  the  850,  700  and  900  mb  meso-scale  analyses.  The  diagnosed 


flow  at  TOm  is  shown  in  i  igure  10.  Tne  super  imposed  isobaric  anaLysis, 
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fhe  12  hour  forecasts  of  winds  at  1600m  and  50m  AGL  are  depicted 
in  figures  11  arid  12,  respectively.  Compared  with  the  independent 
verifying  analysis,  it  is  apparent  that  the  model  showed  a  more  open 
wave  in  the  Great  Lakes  Legion  with  the  warm  frontal  trough  over  Lake 
Huron  and  ^ake  Ontario.  The  model  tended  to  push  the  cold  air  south¬ 
ward  along  the  mountains  more  rapidly  and  eastward  in  Wisconsin  more 
slowly  than  the  verifying  analysis  indicated.  The  flow  forecast  in 
the  vicinity  of  the  high  centered  over  the  four  corners  area  of  Colorado, 
Lew  Mexico,  Arizona  and  Utah  appears  to  be  extremely  accurate. 

The  24  hour  forecast  of  winds  at  i600m  and  50m  are  shown  in 
Figures  13  and  14,  respectively.  It  is  interesting  to  note  that  the 
50m  wind  field  indicates  the  frontal  system  shown  on  the  verifying 
analysis  over  Levada  could  be  nearer  the  Utah  border  (surface  observa¬ 
tions  shown  in  Figure  8  tend  to  support  this  conclusion).  The  fore¬ 
cast  of  the  winds  associated  with  the  leading  edge  of  the  arctic  high 
trailed  the  actual  frontal  position  by  a  full  grid  length  at  the  end 
of  the  24  hour  forecast.  The  lee  side  trough  over  Eastern  Colorado 
with  the  associated  wave  on  the  arctic  front  in  the  I'ebr?'  ka  panhandle 
is  well  represented  by  the  boundary  layer  forecast. 
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•  he  frontal  positions  from  Figures  4,  n  and  8  were  consolidated 
into  a  continuity  chart  shown  in  Figure  li>.  1’his  figure  also  shows 
tne  location  of  grid  row  il  in  relation  to  the  arctic  high  east  of  the 
rocky  Mountains. 

f he  wind  directions  displayed  in  1  igure  16  (  a  time  cross-section 
at  grid  point  near  Medford,  iregon)  are  given  in  relation  to  the  co- 
■rdinate  system  used  in  tne  model,  and  it  is  necessary  to  subtract 
15  fro:.;  tne  direction  to  obtain  wind  direction  in  relation  to  true 
north,  it  is  interesting  to  note  the  development  of  the  wind  max  at 
approximately  300m  above  the  surface  during  tile  4-6  hr  forecast  period, 
associated  witn  frontal  passage  forecast  for  the  area.  The  isotherm 
;  attern  also  Is  quite  indicative  of  the  frontal,  passage  and  shows  the 
relatively  homogenous  air  mass  behind  the  front. 

an  entirely  different  situation  is  depleted  in  Figure  17,  which 
snows  a  time  crcss-section  at  a  grid  point  near  Limon,  Colorado  (sub¬ 
tract  Ip1-'  from  all  wind  direct  ions  to  convert  to  true  directions).  At 
tine  initial  tine,  relatively  warm  air  was  over  tliis  grid  point.  Later, 
the  arctic  air  over  Nebraska  moved  southward  bringing  cold  air  over  the 
gri.d  point ;  a  wave  developed  on  the  arctic  boundary  and  the  cold  air 
receded.  Westerly,  down-slope  winds  then  caused  strong  adabiatic 
warring  which  is  in  evidence  at  tne  end  of  the  forecast  period. 

The  eastward  movement  of  the  arctic  high  with  the  rapid  occlu3sion 
of  the  frontal  system  is  shown  in  'figures  18,  19  and  20.  In  these 
figures,  the  scale  is  greatly  distorted;  the  vertical  coordinate  is 
approximately  two  miles  while  tne  abscissa  is  nearly  3000  miles.  The 
cell  letters  cf  stations  are  only  approximate  locations  and  may  be  either 
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north  or  south  of  the  actual  grid  row. 

-values  computed  in  the  post-processing  phase  of  the  operational 
run  have  proven  very  useful  In  locating  and  predicting  frontal  move¬ 
ments.  nn  interpolated  13-value  field  at  1600m  AGL  is  depicted  in 
rigure  21.  an  analysis  of  the  13-values  computed  for  the  surface  is 
shown  in  l'igure  22.  frontal  positions  together  with  high  and  low  pres¬ 
sure  centers  (from  i  igure  4)  were  added  to  show  correlation  with  the 
surface  analysis. 

The  12  hour  forecast  of  D-values  at  1600m  AGL  (from  the  free  atm  ;s- 
pheric  model)  is  shown  in  !•  igure  23.  The  computed  surface  D-value 
from  the  forecast  temperature  structure  is  shown  in  Figure  24.  The 
boundaries  of  the  arctic  air  are  clearly  defined  on  this  chart.  Fore¬ 
casters  often  do  not  realize  the  great  change  in  D-value  which  can  occur 
vertically  in  the  boundary  layer  region.  For  example,  note  the  arctic 
high  over  Western  Canada  where  the  D-value  at  1600m  AGL  is  approximately 
40  feet,  while  at  the  surface  the  D-value  is  +690  feet. 

Temperature  soundings  for  two  stations  are  shown  in  Figure  25. 

The  sounding  for  Albany  shows  that  in  12  hours  a  slightly  unstable 
layer  between  150  and  600m  became  a  stable  inversion  whereas  at  Omaha 
an  inversion  between  300  and  600m  became  more  unstable  with  the  strong 
advection  of  cold  air.  The  very  strong  gradients  associated  with  an 
arctic  out-break  of  this  type  require  a  very  precise  initial  re¬ 
presentation  of  the  thermal  structure  and  a  very  accurate  forecast 
of  winds  at  all  levels  to  completely  capture  the  forecast  change. 

This  winter  case  study  demonstrates  the  capability  of  the  model 
to  analyze  arid  forecast  detailed  changes  with  considerable  fidelity. 
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It  is  anticipated  that  further  refinement  of  radiation  coefficients 


and  roughness  fields  will  be  made  as  daily  evaluation  of  the  Model  con¬ 


tinues.  Further  development  in  the  area  cf  exchange  coefficients  for 


transport  of  heat  and  moisture  by  eddy  motion,  and  a  more  refined 


thermal  wind  computation  may  further  improve  the  forecasts  from  this 


model. 
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FIG.  16  TIME  CROSS-SECTION  OF  WINOS  AND 
TEMPERATURES  FOR  MEOFORD,  OREGON 
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12  SPONSORING  MILITARY  ACTIVITY 


Air  Force  Global  ’weather  Central 
Cffutt  Air  l'orce  Base,  Nebraska  68113 


i3  abstract  ^  ignited  area  seven  layer  pny sicaa-numericai  inodel  for  the  lower 
tropospnerlc  region  (surface -1600m)  is  described.  The  grid  interval  is  half  that  of 
the  standard  numerical  weather  prediction  grid  used  in  the  hemispheric,  free  atmosphex 
operational  model  at  tne  Air  force  Clobal  Leather  Central  (AICWC).  This  model  is  an 
integral  part  of  the  complete  A] GLG  meso-scaie  (sub-synoptic)  numerical  analysis  and 
prediction  system.  This  model  provides  greater  horizontal  and  vertical  resolution 
in  both  tne  numerical  analyses  and  numerical  forecasts.  It  is  used  to  predict  the 
more  detailed  smaller  scale  atmospheric  perturbations  which  are  important  in 
specifying  sensible  weather  elements. 

important  features  of  this  boundary  layer  model  include;  a  completely  automated 
objective  numerical  analysis  of  input  data:  the  transport  of  heat  and  moisture  by 
three  dimensional  wind  fl< w  (including  terrain  and  frictionaily  induced  vertical 
motions),  intent  heat  exchange  in  water  substance  phase  changes;  and  eddy  flux  of 
iieab  and  water  vapor. 

input  data  ire  conventional  synoptic  surface  and  ut.per  air  reports.  Other 
, redaction  models  provide  horizontal  wind  components  at  the  upper  boundary  and  an 
estimate  of  cloudiness  above  tne  boundary  layer  region,  forecasts  for  the  lower 
boundary  and  surface  iayer  are  empirically  derived.  Despite  some  approximations 
wiilcm  broadly  simplify  the  real  planetary  boundary  layer  processes,  operational 
use  indicates  tne  .model  is  capable  of  producing  detailed  forecasts  but  to  24 
hours.  winter  case  study  is  discussed. 
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